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Abstract 


This  paper  is  concerned  primarily  with  determining  the  mutual  coherence 
function  of  the  field  produced  by  a  plane  quasi-monochromatic  source  in  a  re¬ 
gion  of  variable  refractive  index.  A  scalar  theory  is  used  throughout. 

A  brief  review  of  the  conceptual  background  of  coherence  theory  is  pre¬ 
sented  in  Section  1.  Section  2  contains  an  outline  of  the  mathematical  formal¬ 
ism  of  coherence  theory  and  shows  that  in  a  region  of  variable  refractive  index 
the  mutual  coherence  function  is  propagated  according  to  a  pair  of  inhomoge¬ 
neous  scalar  wave  equations.  In  Section  3  the  pair  of  wave  equations  are  solved 
with  appropriate  Green's  functions.  An  expression  is  derived  for  the  mutual 
coherence  function  of  a  field  produced  by  a  plane  quasi-monochromatic  source. 
In  Section  4  the  case  of  a  statistically  homogeneous  medium  is  treated  and  an 
expression  for  the  ensemble  average  of  the  mutual  coherence  function  is  obtained 
in  terms  of  integrals  of  the  two-point  correlation  function  characterizing  the 
medium. 


ON  THE  MUTUAL  COHERENCE  FUNCTION 
IN  AN  INHOMOGENEOUS  MEDIUM 


1.  Introduction 

The  development  of  coherence  theory  has  been  strongly  influenced  by  re¬ 
search  in  visual  optics.  The  theory  is  concerned  with  the  behavior  of  electro¬ 
magnetic  fields  at  frequencies  so  high  that  measurements  can  be  made  only  by 
averaging  intensities  over  periods  of  time  that  are  long  compared  with  the  times 
involved  for  individual  fluctuations  of  the  fields.  The  fields  are  assumed  to  have 
stationary  time  dependence,  at  least  for  the  intervals  of  the  averaging  periods. 

In  an  advanced  theory  of  optical  behavior,  the  concepts  of  amplitude  and 
phase,  which  are  often  helpful  in  other  branches  of  electromagnetic  theory  and 
in  elementary  optics,  are  no  longer  useful.  The  usefulness  of  these  concepts 
breaks  down  not  only  because  the  high  frequencies  involved  in  optics  make  it  im¬ 
possible  to  measure  the  amplitude  and  phase  of  field  components,  but  even  more 
fundamentally,  light  in  its  usual  form  (the  superposition  of  a  large  number  of 
randomly  timed  statistically  independent  pulses)  is  not  strictly  monochromatic 
but  consists  of  spectra  of  finite  widths.  Indeed,  under  the  conditions  stated  it 
cannot  be  analyzed  even  by  a  Fourier  decomposition  into  strictly  monochromatic 
components  because  only  the  power  carried  by  a  narrow  band  of  wavelengths 
can  be  measured. 

In  the  more  elementary  theory  of  optical  behavior,  the  concepts  of  ampli¬ 
tude  and  phase  presuppose  a  strictly  monochromatic  source  of  illumination  and 
are  not  measurable  quantities  (at  least  in  the  optics  realm).  The  basic  quantities 
in  coherence  theory  (time-averaged  intensities  and  functions  that  express  the 
degree  of  correlation  between  the  vibrations  at  different  points  in  the  field), 
however,  are  measurable  and  do  not  presuppose  the  nature  of  the  field.  Further¬ 
more,  as  limiting  cases,  coherence  theory  yields  not  only  a  strictly  monochro¬ 
matic  theory,  but  also  a  theory  of  incoherent  radiation  (addition  of  intensities) 
and  gives  anaccurate  description  of  the  region  of  partially  coherent  light  between 
these  two  limits. 
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2.  Outline  of  Mathematical  Formalism  of  the  Theory 

The  fundamental  entities  of  coherence  theory  —  the  mutal  coherence  func¬ 
tion  and  the  complex  degree  of  coherence  —  will  be  introduced  in  this  section, 
and  the  wave  equations  for  the  propagation  of  the  mutual  coherence  function  in 
an  inhomogeneous  medium  will  be  derived.  For  the  most  part,  theorems  will 
be  stated  and  the  reader  referred  to  the  literature  for  proofs  and  further  dis¬ 
cussion. 


In  a  source-free,  although  not  necessarily  homogeneous,  medium  the  real 
scalar  function  of  position  and  time  Vr(P,  t)  satisfies  the  scalar  wave  equation 

i  82Vr(P,  t) 


V2Vr(P,  t)  = 


C2(P) 


at 


where  P  is  the  position  vector.  (It  will  occasionally  be  convenient  to  indicate 
position  by  a  subscript  or  to  omit  explicit  spatial  dependence. ) 

The  intensity  of  Vr(t)  averaged  over  an  interval  of  time  of  length  2T  is 
given  by  the  expression 


IP,  t)  ]2dt. 


In  the  following  applications,  T  will  be  extremely  large  in  terms  of  time  units 
of  the  order  of  the  actual  fluctuations  (for  example,  the  mean  period  l/v  where 
v  is  the  mean  frequency  of  the  disturbance) .  It  is  therefore  convenient  to  let 
T— ►<*>  in  expressions  for  the  time-averaged  intensity  and  for  the  correlation 
functions  to  be  introduced  below.  It  is  assumed,  of  course,  that  the 


lim  —  CT  [Vr(P,  t)]2dt 
T-*»  2T  '-'-T 


is  finite. 


A  complex  function  V(P,  t),  the  analytic  signal,  is  associated  with  the  real 

function  Vr(P,  t).  The  advantages  of  choosing  the  analytic  signal  as  a  complex 

o 

representation  of  the  disturbance  have  already  been  discussed  at  length.  The 
analytic  signal  may  be  defined  in  terms  of  Hilbert  transforms  as  follows: 


Let  fr(t)  be  a  real  function  of  time  such  that  its  Hilbert  transform  exists 
(a  sufficient  condition  is  square  integrability);  the  analytic  signal,  f(t)  [associ¬ 
ated  with  fr(t)  j,  is  defined  to  be  f(t)  =  fr(t)  +  if*(t),  where  f*(t)  is  the  Hilbert 
transform  of  fr(t);  that  is. 


At) 


2l 

77 


fr(t') 


t'-t 


dt'sH  [fr(t)  ] . 
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The  integral  sign  j*  denotes  that  the  Cauchy  principal  value  is  taken  at  t'=t. 

Although  Vr(t)  is  in  general  not  square  integrable,  it  is  possible  to  employ 
an  analytic  signal  representation  by  the  stratagem  of  truncating  Vr(t)  at  a  particular 
value,  t=T,  and  then  let  T— »  »  .  Thus,  let 


Vr(P;  T,  t)  = 


Vr(P,  t)  |t|  <  T 
0  |t|  >  T 


and  the  associated  analytic  signal  V(P;T,  t)  be  given  by 


V(P;  T,  t)  =  Vr(P;  T,  t)  +  ^V^P,  t) 


where 


^(P,  t)  =  H[Vr(P;T,  t)]. 

It  should  be  noted  that  ^,V1(P,  t)  is  not  a  truncated  function. 


In  terms  of  these  functions  the  basic  quantities  of  coherence  theory,  the 
mutual  coherence  function,  and  the  complex  degree  of  coherence  may  be  pre¬ 
cisely  defined.  The  mutual  coherence  function  Pg,  T)  =  g^7 )  is  defined 

as  the  complex  cross  correlation  between  the  analytic  signal  representation  of 
the  real  field  at  the  two  points  Pj  and  Pg ;  that  is. 


ri2<T> 


_  lim  f 
'  T-« 


V^T.t  +  r)  V|  <T,  t)  dt. 


The  time  average  and  limiting  proc  ess  is  denoted  by  sharp  brackets: 

r!2(T)  =  <Vt+T)  V2(t)>* 


By  use  of  the  theorem  that  the  cross  correlation  of  the  two  real  functions  is 

'  3 

equal  to  the  cross  correlation  of  their  Hilbert  transforms  in  the  same  order, 

it  can  be  readily  shown  that  the  time-averaged  intensity  Im  at  the  point  Pm  is 

given  by  1/2  rmm(0);  that  is, 

1„  ■  <!vmW>  ■  !/2  W0'-  <»  ■ 1 21  • 

Also,  by  use  of  the  theorem  that  the  convolution  of  two  analytic  signals  is  itself 
an  analytic  signal,  it  can  be  shown  that  is  an  analytic  signal.  From  the 

last  result  it  follows1  that  possesses  a  Fourier  spectrum  that  is  zero 

for  half  the  frequency  range.  Thus 


r12<r>-  r 

OO 
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where 


•r 


ri2(T)  e2*iVTdr  -0P<0. 


It  should  be  noted  here  that  in  formulating  the  solution  to  actual  problems, 
expressions  of  the  form  (tj  +  t)  v2*(t2  +  t)]>  are  usually  obtained.  Under 

the  change  of  variables  t'  =  t  +  tg ,  <  V^t1  +  r)  Vg  ♦(t1)^  is  obtained,  where  r  = 
tg-tf.  The  additional  assumption  of  stationarity  of  V(t)  (that  is,  the  time  aver¬ 
ages  are  independent  of  the  choice  of  time  origin,  or  equivalently  that  the  time 
averages  are  a  function  of  the  difference  in  time  only)  is  necessary  to  equate 
/^2<t)  with  ^Vjftji+t)  Vg  *(t2  +  t)^>. 

A  normalized  form  of  the  mutual  coherence  function,  called  the  complex 
degree  of  coherence  and  denoted  by  y  ^2(t),  *s  verT  useful  in  coherence  theory; 
y  12<r)  *s  defined  to  be 


12 


(  r) 


rl2(T) 


Vrn(0)  r22< o) 


With  the  Schwartz  inequality  it  can  be  shown  that  0<|y12(r)  |  <.1.  The  limits 
characterize  incoherent  radiation  and  coherent  radiation  respectively. 

The  important  result  that  /^(r)  is  propagated  according  to  a  pair  of  wave 
equations  will  now  be  proved.  Specifically: 


V^l'  S2.  *  ' 


c2  (P  ) 

~m 


a2  r(Pv  P2,  t ) 

IP 


-,  (to  =  1,  2).  (2.  1) 


Here  the  Laplacian  acts  on  tne  coordinates  of  the  point  Pm  (m  =  1,  2),  and  tne 
spatial  dependence  of  tne  velocity  of  propagation  is  indicated. 

To  prove  Eq.  (2.  1),  it  is  assumed  that  Vr(P,  t),  hence  the  truncated  function 
Vr(P;  T,  t)  satisfies  the  scalar  wave  equation 


V2Vr(P;  T,  t)  = 


C2(P) 


32  Vr(P;  T,  t) 

at2 


(2.2) 


Also,  as  defined  above 


.V1  (P,  t) 


Vr(P;  T,  t') 
t'-t 


dt'. 


(2.3) 


Both  sides  of  Eq.  (2.  3)  are  operated  on  with  the  Laplacian  V  ,  the  order  of 

4 

operations  is  interchanged,  and  with  the  theorem  that  the  Hilbert  transform  of 
the  derivative  of  a  function  equals  the  derivative  of  the  Hilbert  transform  of  the 
function,  the  following  is  obtained: 


5 


Vrrf  V1  (P,  t)  = 


!  3x2Vi(P.  t) 


1  “  CT  (P)  at* 

Multiplication  of  Eq.  (2.  4)  by  i  and  addition  to  Eq.  (2.  2)  yields 


V2V  (P;  T,  t)  =  . 

CT(P)  atz 

Thus,  the  analytic  signal  itself  satisfies  the  wave  equation. 


r(Pr  p2,t)  = 


—  C  V(P,;  T,  t +r)  V*  <P,;T,  t)  dt.  (2.6) 
2T  J-»  1  z 


Equation  (2.  6)  is  differentiated  with  respect  to  Pp  the  order  of  operations  is 
interchanged,  and  Eq.  (2.  5)  is  substituted,  yielding 


Vi2  r12(r)  4^  -L  £L  VnvT.t  +  rnVcr, 


lim  i  C  1  ^  Vj  (T,  t  +  r ) 

—  T-»  oO  „  9  o 


c2^) 


V2  (T,t)dt 


~2 - 5T  ^12^" 

CT^)  8r2  12 


Similarly, 


^2  ^12^ 


1  3  ri2^> 

c2(p2)  ar2 


and  the  proof  is  complete. 

3.  Mutual  Coherence  Function  in  an  Inhomogeneous  Medium 

The  central  problem  of  this  paper  is  the  determination  of  the  mutual  co¬ 
herence  function  for  a  field  produced  by  an  extended  polychromatic  source  in  a 
region  of  variable  index  of  refraction.  In  the  discussion,  S  is  an  arbitrary  sur¬ 
face  containing  an  extended  polychromatic  source  with  a  known  distribution  of 
mutual  coherence;  P^  and  P2  are  points  in  the  illuminated  field  V;  and  and  S2 
are  points  on  the  surface  S. 

As  shown  above,  the  propagation  of  the  mutual  coherence  function  in  a  source- 
free  but  inhomogeneous  medium  is  governed  by  the  pair  of  wave  equations 


V  2  r  (r)  =  _ - _ 

Vm  'l2(r)  c2(p  } 

— m 


3‘  ri2(r) 


,  m  =  1,  2. 
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It  is  assumed  that  7"^2(t)  is  known  f°r  all  pairs  of  points  and  Sg  on  the  surface 
S. 

Let  ^12^  the  Fourier  transform  of  /"^(t).  Then»  as  stated  above, 
since  /|2  (t)  is  an  analytic  signal,  its  Fourier  spectrum  contains  positive  fre¬ 
quencies  only;  that  is. 


ri2{r)-f  P12(v)e'2Til'Tto  (3.2) 

where 

ri2M  =  ri2(T)e2*iuT  dr.  (3.3) 

Substitution  from  Eq.  (3.  2)  and  Eq.  (3.  1)  and  interchange  of  the  order  of  in¬ 
tegration  and  differentiation  yields 

Jo  [Vm2  +  km2(l;)J  ^12(l,>e"2,rlI/T  d"  =  0  (m  =  l.  2).  (3.4) 

Since  Eq.  (3.4)  holds  for  all  t. 


[V2  +  k2(Pr  V)] 

O 

II 

CO 

M 

A 

(3.  5a) 

[V22  +  k2(P2,  V)\ 

r(Px ,  p2,  v)  =  o. 

(3.  5b) 

Here  k  (P,  v)  =  2ttv/C(P).  Thus,  each  spectral  component  of  V"^2(r  )  satisfies 
the  pair  of  Helmholtz  equations,  Eqs.  (3.  5a  and  b). 

In  Eq.  (3.  5a),  Pg  is  a  fixed  parameter  so  far  as  the  operator  is  concerned. 
In  particular,  Eq.  (3.  5a)  holds  if  Pg  is  a  fixed  point  £>gOn  a  closed  surface  S. 

The  equation  then  becomes 

[V!2  +  k2(Pr  V)  ]  HPj,  S Vv)  -  0.  (3.6) 

The  boundary  condition  for  Eq.  (3.  6)  is  the  known  values  of  r(Sy  S2,u).  Hence, 
the  problem  is  to  solve  the  pair  of  equations 

1*1 2  +  k2(Px,  riinPj,^,!/)  =0  (3.  7a) 

A 

with  £>  2  •  ^  known  on  the  boundary,  and 

[v22  +  k2(P2,  „)j  npv  p2,V)  =  o 


(3.  7b) 
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with  f(P  ,  S  g ,  v)  known  on  the  boundary  as  a  result  of  solving  Eq.  (3.  7a). 

A  formal  solution  to  Eqs.  (3.  7a  and  b)  can  be  obtained  easily  in  terms  of 
the  Green's  function  G  (P,  P')  which  satisfies  the  equation 

[V2  +  k2  (P)]  G(P,  P')  =  -6  <P-P')  (3.8) 

and  which  vanishes  on  the  boundary  S.  In  the  case  that  S  is  plane,  the  Green's 
function  must  be  chosen  not  only  to  vanish  on  the  boundary  plane  but  also  to 
satisfy  the  radiation  condition  at  infinity. 

In  the  same  way  as  if  k  were  a  constant,  the  following  is  obtained: 


"f  I 


r(Si's2') 


h'  Si' 


9Gi<Pi.  sx’) 


8G2(P2,  S  ') 

*n2' 


dSj'  dS2 


(3.9) 


In  Eq.  (3.  9)  and  in  following  equations  primes  are  used  to  indicate  the 
variables  of  integration,  and  corresponding  subscripts  and  primes  attached 
to  the  surface  S  also  serve  this  purpose.  Explicit  dependence  on  the  frequency  v 
has  been  omitted  for  conciseness.  It  will  be  shown  below  that  Gg  =  G^*. 

An  explicit  form  for  G  (P,  P')  can  be  obtained  using  an  iterative  procedure. 

p 

Equation  (3.  8)  is  rewritten  to  represent  k  (P,  v)  as  the  sum  of  a  fixed  mean 

2  ~  p 

value  k  ( v )  and  a  variable  part  with  zero  mean,  k  Me  (P,  i>).  Thus, 

(V2  +  k2)  G(P  ,  P')  =  -6(P-P’)-k2e(P)G(P,  P').  (3.  10) 

(Physically  e(P)  can  be  said  to  be  the  variable  part  of  the  dielectric  constant  of 
the  medium.)  Now  Eq.  (3.  10)  can  be  formally  taken  to  be  an  inhomogeneous 
constant- coefficient  Helmholtz  equation  with  the  right-hand  side  as  the  source 
term.  Acc  -dingly,  a  solution  to  E(p  (3. 10)  can  be  obtained  in  terms  of  the 
Green's  fui,  n  g(P,  F")  which  satisfies  the  constant  coefficient  equation 

(V2  +  k2)  g(P,  P")  = -6(F-P")  (3.11) 

and  vanishes  on  the  boundary  surfaces  Since  g  a'lso  vanishes  on  the  surface, 

G(P,  P')  =  g(P,  P ')  +  k2  j1  e  (P  ")  G(P  ",  P  ')  g(P,  P  ")  dV".  (3  12) 

V" 

Equation  (3. 12)  can  now  be  used  as  the  basis  for  an  iterative  development: 

G(P,  P')  =  g(P,  P')  +  k2  C  e(P")  g(P",  P ')  g(P",  P  )  dV" 

vM 

arc-  (3- 13) 

+  k*  \  \  e(P")  e(P"')  g(P"',  P")  g(P",  P')  g(P",  P)  dV"'dV'!' 

y"  ym 

+  .. . 

Substitution  of  the  iterative  series  Eq.  (3. 13)  into  Eq.  (3.  9)  yields  the  following 
expression  for  the  Fourier  transform  of  the  mutual  coherence  function: 


8 


£(d  _  .  f  fA  agAP  S  ')  3g  (P„,  S') 

-1'  -2  =  t  [/-''• v  ^r~  ' d-2' 

'*{.  If*'- v-v.  p,.  ^  -v  *  *■ 

*  *’  f  .  f f ' -1’’  *■’  .-*.■•  H*"-  W  IV  as,1  as2. 


-2  -1 


+k4I  pgr.v— f  «f r> 

Sj'Sj'  %  v^'  v  "• 

•  -i'"’  gjT,". 

3g. 

,  P,)— i  dV,"'  dV,"  dS,'  dS  ’ 

aijj'  ”*  -1  ~2 

*k*l,  -21’ l  „  I  ‘  <Tl">  f  (f2")  P,)  g2(P2".  P2> 

—2  — 1  V1  v2 

•*i 

*2  <P2"'  V> 

*  <-  — 2  _it f  il  'I  jo  i  jp 

,  —2  -1  d-2  -1 

+  fc«f  f£,=  ,  c  „*i  “VVr  c 

H  1  lAr  2  -2'",‘2-2’" 

.  P2">g,<P2'’. 

?2>  5j’  d"2'"  d-2"  d-f  «l' 

(3.  14) 


The  Green's  function g  (P,  P' )  has  been  determined2  for  the  important  case  in 
which  S  is  a  plane  surface.  The  result  is  given  here: 


g(P,  P')  = 


|P-P'I 


JP-P' 


The  image  of  the  point  P  in  the  plane  S  is  denoted  by  P.  (see  Fig.  1).  In 
Eq.  (3. 14)  the  plus  sign  is  taken  where  a  subscript  1  appears,  and  the  minus 
sign  is  taken  where  a  subscript  2  appears. 


It  should  be  noted  that  an  identical  expression  for  F(Plf  P2 ,  p)  can  be 
obtained  by  a  procedure  that  iterates  for  the  transform  of  the  mutual  coherence 
function  itself  rather  than  for  the  Green's  function  of  the  Helmholtz  equation 
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with  a  variable  coefficient.  This  procedure  puts  Eqs.  (3.  7a  and  b),  the  pair  of 
differential  equations,  in  the  form  of  a  pair  of  equivalent  integral  equations  by- 
use  of  Green's  functions.  Thus, 


A  p  A  agJP-.S') 

r  (Pj.  s2,  v)  =  -  j  r< s1>,  s2,  v )  ds1 

R  i  8*7 1' 


§1 


+  k2  j1  e  (P^inPj',  S2)  g^,  P^)  dVx' 

, 

A  1  pA  9g9(P9,  S’) 

r(p  p  i/)  =  -  j  r(P  s  ’)  -  2— 2  ~2  ds  ■ 


(3. 16a) 


*2 


+  k‘ 


V  ' 


(3.  16b) 


Here,  the  Green's  function  g  (P,  P')  is  identical  to  the  Green's  function  of  Eq.  (3. 11), 

2  ^ 

and  the  representation  k  (P)  =  k  [1+6  (P)]  has  again  been  employed.  The  in¬ 
tegral  equations,  Eqs.  (3.  16a  and  b),  are  then  used  as  the  basis  for  an  interative 
procedure  that  starts  with  the  surface  terms  as  a  zeroth  order  approximation 
and  by  successive  substitutions  yields  the  same  iterative  series  for  P2 ) 

as  obtained  above. 


The  form  of  the  iterative  solution,  Eq.  (3.  14),  obtained  for  the  propagation 

of  the  mutual  coherence  function  in  a  medium  in  which  the  refractive  index 

> 

varies  is  that  of  the  uniform  space  solution 


r< Sj',  s2-) 


st') 

&nx' 


ds2(P2.  s2') 


dSx'  dS2' 


modified  by  a  series  of  correction  terms  that  involve  volume  integrals  of  a 
quantity  associated  with  the  fluctuations  of  the  refractive  index.  When  the 
refractive  index  is  constant,  the  correction  terms  become  zero  and  the  solution 
reduces  to  the  uniform  space  solution. 

The  iterative  solution  obtained  is  that  for  a  single  spectral  component  of 

A 

r  12(t).  To  obtain  2 ( t ) ,  the  iterative  solution  /^2(r)  must  be  substituted 

in  Eq.  (3.  2)  and  the  integration  over  v  performed.  However,  in  the  case  of 
quasi-monochromatic  fields,  it  is  unnecessary  to  perform  the  actual  integration. 
A  quasi-monochromatic  field  is  one  for  which  the  effective  spectral  range,  bu, 
is  small  compared  with  the  mean  frequency,  j7;  that  is,  Civ/v  <  <  1.  For  this  case, 
provided  that  the  time  difference  r  is  small  compared  with  the  coherence  time 
l/Ay,  it  is  known  ’  that  the  mutual  coherence  function  is  of  the  form 
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ri2(T)  ~  ri2  <°>e'2,ril'T.M<<  1/Av  (3-17) 

where  v  is  the  mean  frequency  of  Substituting  from  Eq.  (3,  17)  into 

Eq.  (3.  1)  yields 

[Vm2  +  k2(Pm,D)]  ri2( 0)  =  0  (m  =  1,  2)  (3.  18) 


where 


k  (P 


■m’ 


v)  = 


2irv  . 

C(P  ) 
— m 


Thus,  under  the  quasi-monochromatic  approximation  of  narrow  spectral 
width  and  small  path  differences,  it  has  been  shown  that  /"^(O)  satisfies  the  same 
pair  of  Helmholtz  equations,  Eq.  (3.  5)  ,  as  does  /^(v)  with  v  now  fixed  at  the 
mean  frequency  v.  The  boundary  condition  becomes  ns^,  S2 ,  0)  and  the  iterative 
solution  yields  7^2(0)  which  when  substituted  in  Eq.  (3.  17)  yields  the  quasi- 
monochromatic  solution.  (It  should  be  noted  that  the  varying  propagation  ve¬ 
locity  implies  that  the  time  difference  r  may  be  different  for  two  paths  of  the 
same  geometric  length. ) 

The  Green's  function  formulation  of  the  solution  in  Eq.  (3.  9)  leads  directly 
to  an  important  result  which  we  state  as  a  theorem. 


Theorem 

The  field  produced  by  a  quasi-monochromatic  coherent  source  extended 
over  a  surface  in  contact  with  a  source-free  time- invariant  medium  with  arbitrary 
refractive  index  variation  (in  space)  is  itself  coherent. 


Proof 

2 

The  proof  of  this  theorem  starts  with  the  result  that  a  quasi-monochro¬ 
matic  field  (source)  is  coherent  if,  and  only  if,  the  mutual  intensity  /"^(O)  can 
be  represented  as  the  product  of  a  wave  function  U,  evaluated  at  P^,  with  it’s 
complex  conjugate  U*,  evaluated  at  P2;  that  is 

r(Py  P2,  0)  =  U(PX)  U*  (P2). 

In  the  quasi-monochromatic  approximation  it  is  seen  from  Eq.  (3.  9)  and  the 
discussion  above  that 


-2wivT  C  C  3G  (P  S  >)  3G  (P  S  ') 

T(P.,  P,,r)~  e  2lrlVT  \  \  r(S.',  S',  0)  — i— i i-  —i—i — 2— 

"2  i'i,  -1  “2  81?,'  9V 


2*  h' 


(3. 19) 
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The  assumption  that  the  source  is  coherent  allows  the  factorization 


r<pr  p2,  o)  = 


r  aG-tP.,  S') 

\  U(S')  - 1  dS' 

i,  W 


U*  (S') 


aG2(P2,  s1) 


-dS' 


&n' 


(3.  20) 


where  the  notation  has  been  changed  slightly  to  make  the  relationship  between 
the  bracketed  quantities  more  apparent.  The  proof  will  be  complete  if  it  can  be 
shown  that 


9Gg  SG* 


dT]' 


W 


3  . 


To  prove  this  relationship  the  following  lemma  is  used: 

r*(P1,  p2.  o)  =  ra?2,  p1#  o). 

Substitution  from  Eq.  (3.  20)  into  Eq.  (3.  21),  yields 


C  U*(S') 


3G  *(P..  S') 

1  ■■  ~  dS' 


5 


9G  *(P„,  S') 


■  x  2  -2*  r 

drj' 


I  “S' 


aG,  (P„,  s1) 
_ I  _ 


dS' 


QTi' 


y  u*(s'> 

S' 


0G2(Pr  S') 
377' 


dS' 


(3<  21) 


(3.22) 


V 


p, 


FIG.  1.  Geometry  For  Green's  Function. 
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Since  each  of  the  bracketed  quantities  is  a  function  of  one  point  only,  the  cor¬ 
responding  quantities  can  be  equated,  yielding 

r  9G,*  r  9G„ 

J  U*(S')-^,  (P1,S,)dS'  =  j  U*(.S')  -^-(Pj.^dS1  (3.  23) 

S'  S' 


and  the  desired  relationship  follows  immediately. 

4.  Statistically  Homogeneous  Medium 

An  important  case  is  that  in  which  the  refractive  index  is  a  stationary 
(spatially  homogeneous)  isotropic  random  process.  Assume  that  the  fluctuations, 
e  (P),  satisfy  the  relation 


e(Pi)e(P2)  =  e2(P)  C  (p) . 

Here  e  (PA  e  (P^)  denotes  the  average  of  e  (P^)  e  (P,,)  taken  over  all  pairs  of 

2 

points  P^  and  P 2  a  fixed  distance  p  apart;  €  (P)  is  the  mean  square  deviation 
(of  the  dielectric  constant);  and  C(p)  is  a  correlation  function  that  depends  on 
the  separation  distance  only. 

In  practice,  the  most  that  can  be  expected  from  this  formulation  of  the  pro¬ 
pagation  problem  is  a  prediction  of  effects  'on  the  average.  '  For  example, 
suppose  there  is  a  plane  quasi -monochromatic  light  source  on  a  slab  of 
ground  glass  and  it  is  desired  to  predict  the  distribution  of  the  mutual  coher¬ 
ence  function  on  the  far  side  of  the  slab.  The  solution  will  yield  a  prediction 
of  the  coherence  function  averaged  over  a  large  number  of  different  slabs  of 
glass  with  the  same  statistical  properties.  In  other  words,  the  ensemble  aver¬ 
age  distribution  of  the  mutual  coherence  function  can  be  predicted,  but  the  dis¬ 
tribution  for  a  particular  slab  cannot. 

Suppose,  then,  that  a  series  of  measurements  are  made  to  determine  the 
average  value  of  /"  (P Pg,  r)  which  corresponds  to  a  series  of  independent 
but  statistically  identical  samples  of  a  medium  with  refractive  index  fluctua¬ 
tions.  The  source  distribution  /^(S^,  Jg,  t)  and  the  geometrical  relations 
are  assumed  to  be  identical  for  the  entire  series  of  measurements.  The 
average  cam  be  idealized  by  letting  the  number  of  experiments  N  become  very 
large  and  determining  the  limit  of  the  average  as  N — ►  <» .  The  jth  measure¬ 
ment  is  denoted  by  J" ^  (P^  Pg,  r)  and  the  ensemble  average  of 
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f(P1,  P2»t>  .[rtPj,  P2,  r  )j  is  defined  by 

N 

[n?l'  -2’  T^]  =  N->oo  R  ^  rj  (-l’  -2*  T)‘ 

j  =  1 

With  reference  to  Eq.  (3. 14) ,  recalling  that  only  e  varies  from  one  measure¬ 
ment  to  another,  and  with  the  quasi-monochromatic  approximation,  the  follow¬ 
ing  equation  can  be  written: 


ir<P1.P2,T>1^-2'1"{  J  J  re,-,  s2'.  o)  i 

l  S2'  Sj'  ^1 


*1  ag2 


dT]r 


,dSl'  dS2' 


S  n§,'.52'.o)^  J  [«(p1"»gl(p1",p1) 

-  ‘”'2  V, " 


S2*  S1 


x  88 1  <Pl"»  dv^„  dg^t  dg2, 

8J?!1 


k2I  5"  0>^i  J  I .  (P2"»  g2(P2",  P2) 


V  Sl' 


3>?1  y^n 


ag2(P2">  S  *) 

x  — t  i —  dv2"  dS^'  ds2- 


where  the  ensemble  average  has  been  taken  inside  the  integrals.  Since  the 
samples  of  the  medium  associated  with  the  series  of  experiments  are  assumed 
to  be  independent  and  statistically  homogeneous,  the  ensemble  average  of 
€  (Pj")»  e  (P-^")  €  (Pg").  and  so  on  (formed  with  fixed  points  and  different 
samples  of  the  medium)  can  be  equated  with  the  averages  of  these  same 
quantities  obtained  with  a  particular  sample  and  varying  the  points  (distance 
relationships  are  preserved  when  these  enter  into  consideration).  Hence, 


[((Pm"*!  =  e(P)  =0  ,  (m  =  1,  2)  , 

(  €<Pm")e  <Pm,!')]  =  e(Pn")€(Pm'")  =  e  2(P)  C(p)  ,  (m  =  1,  2)  ,  and  so  on. 
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Thus,  to  second  order  terms  the  result  that  the  ensemble  average  of  the 
mutual  coherence  function  (in  the  quasi-monochromatic  approximation)  is 
given  by 


[HPy  P2,  r)]  =  e 

tk4eV)  y  y 


c  i  c  | 
b2  bl 


X  I  I  C<|P1",-P1"|)g1(P1,"P1")  ij  (Pi”.  P,)— ,~1') 

xr  'I  rr  it.  11  1  1  — 1  — 1  I 


V1"  V  •" 


dT)l 

X  dvn<  dVl",  dSj'  dS2' 

+  k4  e2(P)  f  ^  y  HSj',  S2'.  0)  y  y  C(|  P1-'-P2”| )  g^". 


J2  ** 

X  *2^2",?!' 

+  k47^ 


v  "  V  " 

1  v  2 

„  .  ag,  (P/’.  S')  ag  (p  ",  s  •) 

g9(P0  p0)  — - - 1 - -  - - -  dV  "  dV  "  dS.  -  dS  ■ 

drj  <  377  ’  1  2-1-2 


y  c(|p2...-p  »|,  x 

Og  '  1  U  II  TT  III 


1  V"  V, 


2  2 


ag2  (p2"s2-) 


x  82^2’"'  ?2M)  S2(P2".  P2)  — - ^-^dV,'"  dV„"  dS, '  dS„'. 


2  “”1  u“2 


The  statistics  of  the  medium  enter  into  this  expression  in  the  mean  square  of 

the  refractive  index  fluctuation,  €  2 (P),  and  the  two-point  correlation  function 
C(p)  which  must  be  integrated  over  the  volume  of  the  medium. 
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